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ABSTRACT 

The highly flattened distribution of satellite galaxies in the Milky Way presents a 
number of puzzles. Firstly, its polar alignment stands out from the planar alignments 
commonly found in other galaxies. Secondly, recent proper motion measurements re- 
veal that the orbital angular momentum of at least 3, and possibly as many as 8, of 
the Milky Way's satellites point (within 30 degrees) along the axis of their flattened 
configuration, suggesting some form of coherent motion. In this paper we use a high 
resolution cosmological simulation to investigate whether this pattern conflicts with 
the expectations of the cold dark matter model of structure formation. We find that 
this seemingly unlikely set up occurs often: approximately 35% of the time we find 
systems in which the angular momentum of 3 individual satellites point along, or close 
to, the short axis of the satellite distribution. In addition, in 30% of the systems we 
find that the net angular momentum of the 6 best aligned satellites lies within 35 
degrees of the short axis of the satellite distribution, as observed for the Milky Way. 



1 INTRODUCTION 



The system of bright satellites orbiting in the halo of the 
Milky Way has a number of properties that seem unusual. 
Most of these satellites define a tight plane in the sky whose 
normal is roughly per pendicular to the Galaxy's spin angu- 
lar momentum vector llLvnden -Bcll 19761; iKunkel fc Demersl 
1 19761 ; lLvnden-Beiilll982l). In addition, thi s so-called disc- 
of-satellites l|Kroupa. Thies fc Boibl l2005f) appears to be 
well aligned with the orbital angular momentum of indi- 
vidual satellites, indicating some form of coherent motion 
jMetz. Kroupa. fe Libeskindll2008h . 

The observed, highly anisotropic set-up is not an obvi- 
ous outcome of the ACDM model in which the Universe is as- 
sumed to be composed primarily of cold dark matter (CDM) 
and vacuum energy associated with a cosmological constant 
(A). In this paradigm, structures form hierarchically via the 
merging of small subgalactic units formed out of the un- 
stable overdense peaks of the initial matter density field. 
Satellite galaxies are identified with the substructures that 
survive the violent gravitational processes associated with 



hierarchical growth. High resolution N-body simulations of 
galaxy formation have shown that small bound subhalos can 
survive these merging processes and that many remain as 
distinct substructures embedded in a present day galactic 
dark matter halo (e.g. Klypin et al.lll999l; Moore et al.l|l999l ; 
iDiemand, Kuhlen fc Madaull2007l ; ISpringel et alj|2008r >. 



Before the flood of discoveries of new M ilky Way satel- 
lites i n the Sloan Digital Sk y Survey (SPSS ) l|Willman et al] 
120051 ; iBelokurovl l2006al lbl IZuckeif l2006al lbh. it had been 
known for at least 3 decades that rather than being isotrop- 
ically distributed in space, the 11 most luminous 'classical' 
satellite galaxies of th e Milky Way are all aligned along a 
great circle on the sky dLvnden-Bel l 1976; Kunke Qg Demersl 
Il976| y iMetz. Kroupa. fc Jerienl (|2009h have argued that the 
newly discovered SDSS satellites als o lie on or close to this 
great circle. This 'galactic pancake' l|Libeskind et al.ll2005l ) 
of satellites was found to be aligned roug hly perpendicular t o 
the plane of the disc of the Milky Way (|Lvnden-Belllll982l ). 
This particular set-up may be related to the "Holmberg ef- 
fect," the tendency for satellites of extragalactic systems 
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to be preferentially d istributed alone 
A number of studies (Irlolmberg 1969 



A number ot studies (liolmbcrg 19b9; /.aritsky ■ 
ISales fc Lamba3l2004l : lBrainerdl2005l : lYang et all l 



a specific direction . 
' Zaxitskv et~alll997t 

~~m 



fang et alll2005D have 
detected some alignment of satellite galaxy systems rela- 
tive to t heir hosts, but d isagr ee on the nature of its ori- 
entation. iHolmberd ljl969l l and IZaritskv et~aH (jl997h found 
"Milky Way" type alignments (i.e. polar, with satellites 
avoiding the equatorial regions), while planar alignments 
(i.e. with satellites avoiding the polar regions) have been 
obser ved in the 2-degree- fiel d Galaxy Redshif t Surv ey (2dF- 
GRS, IColless et all200ll) bvlSales fc Lambasl d2004h and, in 
the SPSS, bv lBrainerdl (|2005h and lYang et all (120051). 

iKroupa. Thies fc Boilvl (|2005h argued that the 
anisotropic nature of the distribution of the Milky 
Way's satellite population contradicts the expectations 
of the ACDM model since the simple hypothesis that 
satellite galaxies are a random subsample of the entire 
subhalo population would imply near isotropy. However, 
recent numerical st udies of the formation of satellites in the 
ACDM cosmology (fObeskind et al.ll2005l . 120071 ; iKang et al.1 
120051 ; IZentner et al l l2005h have found that there is not a 



one-to-one match between the brightest satellites and the 
most massive subhalos and indeed that the satellite galaxies 
are not a random subsample of the background subhalo 
population. Instead, the most luminous satellites populate 
the subset of subhalos which were the most massive before 
they were accreted. These studies also suggested that 
the tendency for satellites to align themselves in a great 
pancake may be a generic result of infall along filaments 
that form in the hierarc hical buildup of a galactic halo 
or of accretion in groups (|Li fc Helmill2008l ). This kind of 
coherent motion should perhaps also be reflected in the 
kinematics of the satellite galaxies. 

Proper motions for some of the Milky Way's satel- 
lite galaxies (at the level of milliarcsec /yr) have now 
been measured using both ground- (|Scholz fc Irwin] 



1994 llbataet al.1 Il997l; iDinescu et al.l |2004 120051 : 



Schweitzer et al.l 1 19951. 1 19971) and space-based te le scopes 



(i.e. HST iKallivavalil. van der Marel fc Alcockl 
Kallivavalil et all 120061 ; iPiatek et al.1 12001 120051. 



20071 ). Using these data, iMetz, Kroupa, fc Libeskind! (1200 



2006; 



calculated the orbital angular momentum of 6 satellites. 
Remarkably, they found the angular momentum vectors 
of the satellites to be well aligned, generally pointing in 
the same direction, within 30 degrees, o n the sky. This 
effect was n ot seen in the simulations of (Libeskin d et al.l 
120051 . 120071 ) which were shown to be consistent with a 
random orientation for the satellites' angular momenta. 
This conclusion, however, was based on a very small sample 
of simulations of Milky- Way systems. 

In this study we use a high resolution simulation of a 
large cosmological volume to explore the alignment of the 
positions and angular momenta of satellites in dark matter 
halos. With a much large sample than those previously con- 



Table 1. Parameters of the hMS N-body simulation. (1) Box size 
(in h _1 Mpc); (2) number of particles; (3) dark matter particle 
mass (in Mq); (4) mean matter density; (5) vacuum energy den- 
sity; (6) Hubble constant (in units of 100 kms _1 Mpc~ 1 ); and (7) 
normalization of the matter power spectrum. 



'"■cln 



100 7.29x10 s 1.3x10 s 0.25 0.75 0.73 0.9 



sidered, we are able to investigate these properties with a 
high degree of statistical significance. 



2 METHODS 

2.1 N-body simulation and semi-analytic model 

The parameters of the hMS ACDM N-body simulation that 
we analyze (listed in Table 1) are the same as the Millen- 
nium Simulation l|Springel et alll2005l ) except that the box 
size is smaller a nd the mass reso lution approximately ten 
times better (see lGao et al.l [20081 . for details). Th e simula- 
tion was performed using the L-Gadget2 code (|Springell 
2005) on 128 processors at the Cosmology Machine of 
Durham University's Institute for Computational Cosmol- 
ogy. The simulation was populated w i th 'ga laxies' using the 
semi-analytical model of lBower et al.l ((2006) applied to halo 
merger trees extracted from the simulation. 

Large high-resolution simulations such as ours resolve a 
wealth of gravitationally bound substructures inside larger 
objects. When a small halo is accreted by a larger one, tidal 
stripping causes the substructure to lose mass, while dynam- 
ical friction against the parent halo's background material 
causes its orbit to spiral in towards the centre. In order to 
identify and follow these accreted substr uctures we use the 
algorithm Subfind |Springel et al.ll200ll ). which uses an ex- 
cursion set approach recursively to identify local maxima in 
the dark matter density field. The extent of a substructure is 
then determined by locating saddle points in the dark matter 
density field. The gravitational binding energy is calculated 
and unbound particles are rejected. Finally a 20 particle 
lower limit is imposed on the substructures. 

We identify subhalos at all redshifts in the simulation 
and use them to build halo merger trees w hich describe 
the h ierarchical build up of structures (see lHarker et al.l 
2006, for a full description). The trees provide the back- 
bone for modelling the evolution of the baryons. We use the 
Durham semi-analytic galaxy formation model, Galfo rm 
l|Cole et al.ll2000l; iBenson et al.ll2002l ; iBaugh et aT1l2005l ). as 
extended bv lBower et al.l l|2006h which, in brief, includes the 
following physical processes: (i) shock-heating and virializa- 
tion of halo gas; (ii) radiative cooling of hot halo gas onto 
a galactic disc; (iii) star formation from cool gas in the disc 
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and during starbursts in the bulge; (iv) the evolution of stel- 
lar populations; (v) photoionization and its effect on the 
thermal properties of the intergalactic medium; (vi) reheat- 
ing and expulsion of cooled gas due to feedback processes 
associated with supernovae, stellar winds and AGN; (vii) 
chemical evolution of gas and stars; (viii) reddening due to 
dust absorption; (ix) galaxy mergers; (x) formation of black 
holes in the bulges of galaxies via the accretion of gas during 
mergers and disc instabilities; (xi) the evolution of the size 
of bulges and discs; 

2.2 The satellite sample 

The semi-analytic model produces a galaxy catalogue at 
each simulation output. Resolution tests, in which we in- 
creased the minimum halo mass used in the merger trees, 
indicate that our catalogues are complete for satellite galax- 
ies with K-band magnitude brighter than -8. In the Milky 
Way, this magnitude limit corresponds to the 11th brightes t 
satellite, (Draco, which has M v = -8.8; see lMateolll998h . 
Although the number of known dwarf galaxies orbiting in 
the Milky way halo has increased by at least 15 in the past 
3 years, as new local group dwarf satellites have been dis- 
covered in the SDSS, all of the newly discovered galaxies are 
faint and fall below our magnitude limit. For this reason, we 
do not consider them in this paper. (Note that the newly 
discovered SDSS satellites are spatially biased because of 
the limited sky coverage of the SDSS so it is unclear how 
their properties compare with those of the classical satel- 
lites.) Thus, when comparing with observations of the Milky 
Way satellite system we only consider the 11 classical satel- 
lite galaxies whose luminosity is greater than our magnitude 
limit. 

In order to identify central galaxies that resemble the 
Milky Way, we apply a halo mass cut which excludes galaxies 
residing in haloes of mass < 2 x 10 11 M Q and > 2 x 10 12 M . 
To refine further the search for Milky Way-type galaxies, 
we also apply a luminosity cut to eliminate central galaxies 
fainter than —20 in the V band. We consider only systems 
that contain 11 or more satellites within the halo virial ra- 
dius, r v ir, defined as the distance from halo centre within 
which the mean interior density is 200 times the critical 
density. 

In order to prevent contamination of an alignment sig- 
nal by nearby central galaxies, we impose the isolation cri- 
teria shown schematically in Fig. [1] We eliminate for the 
galaxy catalogue any two central galaxies A and B that lie 
within each other's virial radius (Case I). If galaxy B lies 
within the virial radius of another galaxy A, we eliminate 
both from our sample if B has a stellar mass greater than 
that of A (Case II), and turn B into a satellite of A if it is 
less massive than A (Case III). In practice, cases I, II and 
III are rare with only 29, 63 and 102 galaxies falling into 
these categories respectively. 



Case I 



Case II 



Case III 




Figure 1. Isolation criteria. In case I (top), galaxy A is within 
galaxy B's virial radius and galaxy B is within galaxy A's virial 
radius; in this case both galaxies are eliminated from the sample. 
In case II (middle), galaxy B lies within galaxy A's virial radius 
(but not vice-versa) and galaxy B is more massive than galaxy A; 
in this case both galaxies are eliminated from the catalogue. In 
case III (bottom), galaxy B is within galaxy A's virial radius and 
is less massive than galaxy A; in this case galaxy B is regarded 
as a satellite of galaxy A. 



When applied to the galaxy list produced by the semi- 
analytic model, the halo mass cut leaves us with 30 946 sys- 
tems. The central galaxy magnitude cut reduces our sample 
to 3 230 halos, while the isolation criteria leaves us 3 201 sys- 
tems. This sample is further reduced to 436 galaxies which 
posses at least 11 satellites. The median halo masses of each 
sample after these 4 cuts are 10 1132 , 10 1182 , 10 11 ' 82 and 
lO n ' 96 M0 respectively. Thus the median halo mass is only 
increased by 38% by the requirement that each system host 
at least 11 satellites. 



3 SATELLITE GALAXY ALIGNMENTS 

We first consider the alignment between satellite galaxy sys- 
tems and the shape of the parent halo and then investigate 
the alignment of the satellites' orbital angular momentum 
vectors. 
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3.1 Orientation of the halo and satellite positions 



Once the isolation criteria, and mass and luminosity cuts 
have been applied we obtain a list of central galaxies. We 
wish to examine the flattening of both their haloes and satel- 
lite systems. To do so, we first define the inertia, or second 
moments tensor, 



Ij,k = > Xi 



(i) 



where we sum either over all the particles or all the satellites 
in the halo. We diagonalize this tensor and obtain the princi- 
pal axes of the system, defin ed conventionally such that a > 



b > c. Previous studies (e.g . Bullock. Kravstov fc Weinberg] 
l200d : lLibeskind et alJl2005Ll2007l ) have shown that although 
haloes tend to be mildly aspherical, their satellite popula- 
tions exhibit a much more pronounced flattening. In order 
to verify this with our sample, we calculate the eigenvalues 
of the inertia tensor for the dark matter particles in each 
halo and for the systems of 11 brightest galaxies within r v i r . 

The flattening of the halos in our sample is evident in 
Fig. [2] where we plot c/a versus b/a. From the definition of 
a, b and c, no points may lie in the triangular upper half 
of this plot. The mild asphericity of our sample is reflected 
in the tendency for halos to occupy the upper right area 
of Fig. [2] The median values of these halo axial ratios are 

< b/a >= 0.87 and < c/a >= 0.77. 

The mild asphericity of dark matter haloes may be con- 
trasted with the much more flattened configurations defined 
by their satellite systems, as shown in Fig. [3] This sam- 
ple shows a much broader distribution of axial ratios, and 
contains many more examples of very triaxial systems than 
the halo population. The median values of the axial ra- 
tios are significantly smaller than for the halo sample, with 

< b/a >= 0.69 and < c/a >= 0.43. 

The nature of the satellite's aspherical distribution is 
related to the fact that satellites do not simply trace the 
background dark matter, but occupy biased positions in the 
halo closer to the principal axes. We assess the significance 
of this effect by conducting a "randomization" test of the 
null hypothesis that satellites are randomly drawn from the 
dark matter distribution. For each satellite in a halo, we 
randomly select a dark matter particle at the same radius. 
We then proceed to calculate the principal axes for both the 
satellite systems and those consisting of the random particles 
replacing each of the 11 brightest satellites. 

We use these two sets of principal axes to project the 
location of each satellite and its random counterpart onto a 
Hammer- AitofT projection of the whole sky, shown in Fig. [4] 
The upper panels show the location on the sky of the model 
satellites. For each central galaxy, we plot the positions of 
its satellites rotated onto the reference frame defined by the 
disc-of-satellites of the 11 brightest satellites. By definition, 
the satellites will concentrate along the long axis, a feature 
readily visible in the upper panels of Fig. [4] However, the 



O 




Figure 2. The ratio of short-to-long axis versus the ratio of 
intermediate-to-long axis for our halo sample. Each halo is shown 
as a black point while contours represent the density of points 
on an arbitrary scale, with red being the highest density and 
white the lowest. The large plus sign marks the median value of 
(b/a, c/a) and its arms represent the lcr extent of haloes in either 
direction. The point marked "MW" indicates the location on such 
a plot of the Milky Way's disc-of-satellites 



U 




0.0 0.2 0.4 0.6 



Figure 3. Same as Fig. [2] but for our sample of halo satellites. 
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Figure 4. Upper panel: the location on the sky of all satellites in our halo sample, rotated into the frame denned by the satellite system's 
principal axes. On the left we plot the position of each satellite, while on the right we show a false-colour map of satellite galaxy "density" 
(defined using the angular distance to the 10th nearest neighbor). High density regions are indicated by the purple/red shading, while 
less dense regions are indicated by the white shading. The red, blue, and black circles define areas within 15 degrees of the a, b and, c 
axes. Bottom panel: as above but for "randomized" satellites. 



total fraction of satellites that lie close to the long axis is 
determined by the strength of the intrinsic flattening. 

We can compare the upper panels of Fig. Uto the lower 
panels which show the location on the sky of the random 
particles that replaced the brightest satellites projected onto 
the reference frame defined by those replacing the 11 bright- 
est in each halo. These particles also tend to be aligned with 
the long axis of their distribution (by construction), but to a 
lesser extent than the model satellites. The pseudo- disc-of- 
satellites defined by the randomly selected dark matter par- 
ticles is visibly less aspherical than the true disc-of-satellites 
defined by the 11 brightest satellites in each halo. 

We can quantify this behavior by plotting a histogram 
of the angle between the position vector of each satellite 
(^sat) and the short axis (for example) of the ellipsoid de- 
fined by the brightest 11 satellites (c sa t). Fig.[S^a) shows that 
the model satellites (black lines) exhibit a much stronger 
perpendicular alignment with the short axis of their distri- 
bution than the corresponding set of random dark matter 
particles (red lines). A Kolmogorov- Smirnoff test rejects the 



hypothesis that these two distributions are drawn from the 
same parent distribution or from a uniform distribution at a 
very high confidence level (formal KS probabilities ~ 10 -37 ). 
Around 15% more model satellites lie further than 65 de- 
grees from c sat than expected from the associated dark mat- 
ter particles. We conclude that satellite positions are not a 
random sample of the halo and that satellite systems are 
significantly flatter than their parent halos. 

It is straightforward to check the degree of alignment 
of the disc-of-satellites defined by the 11 brightest galaxies 
with their host halo. Fig. [5jb) shows the cosine of the angle 
between the short axis of the disc-of-satellites (cdos) and 
the short axis of the host halo (chaio)- The model discs-of- 
satellites (black line) show a clear tendency to be aligned 
with the main plane of the halo. Indeed, nearly 40% of all 
satellites systems fall in the last bin of the histogram. How- 
ever, a non-negligible fraction of systems are not aligned, 
in the sense that their minor axes are nearly perpendicu- 
lar to the halo's minor axis. For example, approximately 
10% of the model discs-of-satellites are mis-aligned by more 
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Figure 5. Top panel: (a) The distribution of the cosine of the 
angle between the short axis of the each halo's satellite system 
ellipsoid and the position of each satellite galaxy (black curve). 
The red line shows the analogous distribution for the "random- 
ized" sets of satellite galaxies. Bottom panel: (b) The distribution 
of the angle between the short axis of each halo's dark matter dis- 
tribution and the short axis of the disc- of- satellites composed of 
the brightest 11 satellite galaxies. 



than 75°. To assess the significance of this trend, we define 
randomized discs-of-satellites using the set of dark matter 
particles replacing the brightest 11 model satellites. Again, 
a Kolmogorov-Smirnoff test rejects the hypothesis that the 
actual and randomized samples are drawn from the same 
parent distribution, or that either is drawn from a uniform 
distribution, at a very high confidence level. 



A nu mber 
2001 120051. 



of re cent 
l200rj, [2007 



3.2 Orbital angular momentum of satellites 

studies (e.g. iPiatek et al.l 
Yfl, and other references in 
Metz. Kroupa. fc Libeskind! [2008) have obtained upper 
limits on the proper motion of eight of the Milky Way 
(MW) dwarfs. A detailed study of their angular momenta 
reveal two particularly interesting aspects of their orbital 
motion. The first is that the angular momentum vectors of 
at least 6 of these satellites point to within 30 degrees of 
the mean angular momentum vector of the set for which 
proper motions are available. The second is that the orbits 
of at least 5 satellites fall within 30 degrees of the normal 
of the disc-of-satellites defined by the 11 brightest dwarfs. 



iMetz. Kroupa. fc Lib eskind (2008) tested for these two 
properties in 6 simulations of MW-sized galaxy halos and 
found that in none of them did the satellites (identified in 
a similar manner as in this work) exhibit the same kind 
of coherent motion inferred for a subset of MW satellites. 
None of the 6 simulated haloes in that work had 5 satellites 
with with angular momentum vectors within 30 degrees of 
the normal to the fitted disc-of-satellites. Just 2 of the halos 
had 3 satellites whose angular momentum vectors fell within 
this region. In addition, the mean angular direction of the 
orbital poles of the simulated satellites did not lie close to 
the normal of the disc-of-satellites, in contrast to the MW. 

Metz et al. concluded from this study that the dynam- 
ics of at least 5 "rotationally supported" MW satellites were 
(at least) hard to reproduce in CDM simulations. They sug- 
gested that the mismatch between the simulations and the 
observations could, in theory, be explained if the satellites 
of the MW had a different origin from those in the sim- 
ulations, for example, if they were "tidal dwarfs", formed 
in a tidal stream gen erated in a recent major merger (e.g. 
IMetz fc Kroupall2007r i. 

Our sample of MW type satellite systems is well suited 
to testing whether or not the orbital dynamics of satellites 
in ACDM simulations are consistent with the MW data. 
In Fig. [5] we show the location on the sky of the orbital 
angular momentum of each satellite (upper panels) and its 
randomized counterpart (lower panels) in our halo sample. 
We rotate each point on the sky into the frame defined by 
the principal axes of the "disc-of-satellites" defined either 
by the model satellites or by their randomized counterparts. 
There is an ambiguity in choosing which direction of the c sa t 
axis to take to define the short axis. We choose the c sa t axis 
to point in the direction which is closest to the mean angular 
momentum vector of each system. Therefore, the satellites 
that fall in the Southern hemispheres of Fig. [6] are counter- 
rotating with respect to the mean angular momentum of the 
disc-of-satellites. 

Two features are readily visible in this figure. Firstly, 
the directions of the angular momentum of the randomized 
dark matter particles replacing the model satellites are much 
more uniform than those of the actual model satellites. Sec- 
ondly (by construction), in both samples the distributions 
show an excess around the "North" pole, which is much 
more pronounced for the model satellites than for their ran- 
domized counterparts. On closer examination we note that, 
as well as the higher density around the North pole, there 
is a tendency for the angular momenta to point close to the 
b sa t — S sat plane. This effect may be produced by the radi- 
ally biased orbits of the satellites in our sample. If a satellite 
lies along or near a sa t , and if it is also moving close to this di- 
rection, its angular momentum vector will point somewhere 
on the b sat - 6 sat plane. 

Satellites have a much stronger tendency to display co- 
herent orbital angular momenta then the background dark 
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Figure 6. Upper panel: The location on the sky of the orbital angular momentum vector for all satellites in our halo sample, rotated 
into the principal axes of the disc-of-satellites. In these projections we choose from the two possible opposing directions for the c-axis the 
one which is closest to the mean angular momentum vector of the satellite sample. On the left we show the location of each satellite's 
angular momentum vector, while on the right we show a contour map of satellite angular momentum "density" (defined by the angular 
distance to the 4th nearest neighbor). High density regions are indicated by the red/purple shading, while less dense regions are indicated 
by the white shading. The red, blue, and black circles define areas within 15 degrees of the a, b and, c axis. Bottom panel: The same as 
above, but for randomized satellites. 



matter. In Fig. [7] we show the distribution of the cosine of 
the angle between each satellite's orbital angular momen- 
tum vector, j sa t, and c sat . This plot quantifies what is vi- 
sually apparent in Fig. [5] an above random alignment be- 
tween the orbital angular momentum of individual satellites 
and the short axis of the ellipsoid to which it belongs. The 
Kolmogorov-Smirnoff probability that the model satellites 
and the random particles are drawn from the same popula- 
tion is negligible. 

Thus, our model satellites display the same kind of ap- 
parent coherent motion and spatial flattening as the real 
MW satellite system. Yet, the nature of this coherent mo- 
tion is unclear. Metz et al. suggested that in the MW, the 
disc-of-satellites is supported purely by rotation. In this sce- 
nario, the majority of satellites move on circular orbits such 
that the disc-of-satellites itself is supported by the rotation 
of the galaxies orbiting in it. In order to test this hypothesis 
we carry out the following calculation. As in Fig. [6] we define 



the z-axis to be parallel to the short (c) axis of the satellite 
distribution and of the two opposing directions this defines, 
we select the one closest to the direction of the mean angular 
momentum of the satellites. For each galaxy we calculate its 
angular momentum in this 'z'- direction and call this quan- 
tity Jz(E). We then estimate the angular momentum of a 
satellite with the same orbital energy, but on a purely cir- 
cular orbit, J ZC (E). If these two quantities are similar (i.e. 
Jz(E)/J zc (E) ~ 1), then the satellite's orbit must be ap- 
proximately circular and rotationally supported. Note that 
since a circular orbit maximizes the angular momentum, the 
ratio Jz(E)/J zc (E) must always be < 1. 

For each halo, we estimate its concentration assuming 
a Navarro, Frenk & White (1996; NFW) profile. Once we 
k now the halo concentr ation, we use eqns (2.18) and (2.19) 
of ILacev fc Cole] (|l996T ). together with the satellite's radial 
position in the halo, to calculate its energy. We then search 
for the parameters (i.e. the distance and velocity) of a circu- 
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Figure 7. The distribution of the cosine of the angle between j sa t , 
the satellite galaxy orbital angular momentum, and £ sa t (black 
lines). The red lines show results for the control sample of ran- 
domized dark matter particles. 



lar orbit within an NFW potential with the same energy. In 
this way, we build the ratio Jz(E) / J ZC {E) whose distribution 
we show in Fig. [8] 

There are two salient points visible in Fig. [8] The first 
is that the distributions are not peaked at « 1, but rather 
at » 0, implying that J Z (E) <C J ZC (E). The majority of the 
satellites and their randomized counterparts are therefore 
not on circular orbits around the z-axis. It is interesting to 
note that the control sample exhibits a more pronounced 
symmetry about than the satellite sample, implying simi- 
lar numbers of counter- and co-rotating particles. The satel- 
lite sample, on the other hand, is somewhat skewed towards 
positive values of J Z (E) / J ZC (E) , implying that more satel- 
lites are co-rotating. 

We now quantify how common the configuration of the 
satellite system of the MW is in our ACDM simulation. Of 
the 11 classical satellites in the MW, 8 (Sagittarius, Draco, 
Fornax, Ursa Minor, Carina, Sculptor the SMC and the 
LMC) have proper motion estimates which allow us to esti- 
mate the direction of their orbital angular momentum (e.g. 
see Table 2). Of these eight, the LMC, the SMC, and For- 
nax have angular momentum pointing well within 30 degrees 
of the direction to the normal of the plane fitted to the 11 
classical satellites. The angular momentum vectors of an ad- 
ditional 2, Ursa Minor and Carina, have error bars large 
enough that they are compatible with a < 30 degrees mis- 



0.12 I 




-1 -0.5 0.5 1 

J Z (E)/J Z0 (E) 

Figure 8. Histogram of the ratio of the z- component of the an- 
gular momentum of each model satellite, J Z (E), to the angular 
momentum of a satellite on a circular orbit with the same energy, 
Jzc(E) (black lines). If the disc- of- satellites were rotationally sup- 
ported, all orbits would be circular and the distribution should be 
strongly peaked around ~ 1. The red lines show the the histogram 
for the control sample of matched dark matter particles. 



alignment with the pole. Additionally, there are 3 satellites 
without proper motion estimates, implying a minimum of 3 
out of 11 (~27%) and a maximum of 8 out of 11 (~73%) 
satellites with angular momenta within 30 degrees of the 
perpendicular to the plane of the disc- of -satellites. 

To assess how common this form of alignment is, we 
perform the following calculation: for each simulated halo 
with more than 11 satellites, we fit the standard plane to 
the 11 brightest, as described in Sec. 13.11 We then calculate 
how often the orbital angular momentum of N randomly 
selected satellites falls within 15, 30 or 45 degrees of c aa t, 
where the direction of the c sa t axis is chosen, as before, to 
be the one that is closest to the mean angular momentum 
vector. The results are shown in Fig. [5] 

We find, for example, that the probability of finding at 
least 3 satellites within 30 degrees of c sat is around 35%, 
a similar number to that found in the simulations of Metz 
et al. If n — 5, this value drops to ~ 5%. The probability 
varies rapidly with satellite number and the chance of find- 
ing 9 satellites within the specified angle (the total number 
of classical satellites - save Sagittarius and Sculptor - that 
could possibly be aligned) drops to 1%. The number of satel- 
lites with angular momentum aligned with the normal to the 
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Table 2. Census of the positions and angular momenta of the 
MW's satellite population. Column 1 shows the satellite's name. 
Column 2 gives the galacto-centric distance of each satellite. Col- 
umn 3 gives the angular separation between the normal of the 
plane fitted to the Milky Way's disc- of- satellites and the angular 
momentum vector of each satellite. Column 4 gives angular sep- 
aration between the direction of the mean angular momentum of 
the 6 satellites moving in the most ordered fashion, </ mean , to the 
angular momentum vector of each satellite. Columns 5 and 6 indi- 
cates whether the uncertainty in the measurements of the satellite 
galaxy's angular momentum are consistent with placing the pole 
of the angular momentum vector within 30 degrees of the normal 
of the disc-of-satellite's or within 35 degrees of J^, . Fo r ob- 
servational references sec iMetz. Kroupa. fc Libeskindl ([200^ 1 and 
references therein. 



Name 


d GC 


A0 csat 


A0j 


within 


within 




(kpc) 


(deg) 


(deg) 


"DoS 


</mcan 


LMC 


50 


13 


7 


YES 


YES 


SMC 


57 


21 


15 


YES 


YES 


Fornax 


140 


23 


25 


YES 


YES 


Ursa Minor 


68 


38 


18 


YES 


YES 


Carina 


103 


49 


61 


YES 


YES 


Draco 


82 


69 


55 


NO 


YES 


Sagitarius 


16 


118 


98 


NO 


NO 


Sculptor 


79 


136 


156 


NO 


NO 


Sextans 


86 






unknown unknown 


Leo I 


205 






unknown unknown 


Leo II 


250 






unknown unknown 



disc-of-satellites also depends strongly upon cone angle. For 
example, we find 3 satellites within 45 degrees around 70% 
of t he time but only 5% of the t i me wit hin 15 degrees. 

IMetz. Kroupa. fe Libeskindl (|2008l ) also remarked on 
the coherence of the angular momenta of the Milky Way's 
satellites. They chose to omit 2 satellites (Sculptor and 
Sagitarius) as their angular momentum vectors point far 
from the mean and are clearly not members of a "common 
stream". They then computed the mean angular momen- 
tum vector and its error for the remaining 6 satellites that 
have measured proper motions. Relative to this mean angu- 
lar momentum direction, they found that a minimum of 4 
(the LMC, SMC, Fornax and Ursa Minor) out of 11 (~ 36%), 
and a maximum of 9 out of 11 (~ 82%), satellites fall within 
the 35 degree error circle. 

We perform a similar calculation to determine the prob- 
ability of such an alignment. In order to mimic their proce- 
dure, for each halo we first calculate the net angular momen- 
tum of a randomly selected 8 satellites from the brightest 11. 
Then, after omitting the 2 satellites whose angular momenta 
are furthest from this mean, we compute the mean angular 
momentum vector of the remaining 6. We refer to this quan - 
tity as J mcan . Although IMetz. Kroupa. fc Libeskindl (|2008l ) 
did not use this exact prescription, it is equivalent to theirs 
and it results in the same selection for the Milky Way. We 




2 4 6 8 



satellites 

Figure 9. The fraction of systems in which the angular momen- 
tum of at least N satellites drawn from the brightest 11 lies within 
15 (black), 30 (red), and 45 (blue) degrees of c sat (solid lines) and 
of J^eaxi (dashed lines). 



then calculate how often N randomly selected satellites fall 
within 15, 30 and 45 degrees of the direction of J„ cm . Our 
results are shown as the dotted lines in Fig. [9] 

We find that the arrangement in the Milky Way, where 
4 satellites have angular momentum lying within 30 degrees 
of ./mean, occurs reasonably frequently (~ 10% of the time). 
However, we find that the probability of finding the max- 
imum possible number of satellites (9) with angular mo- 
mentum vectors pointing within the even larger angle of 
45 degrees from J moan is negligible. If future proper motion 
measurements for the remaining 3 satellites (Sextans, Leo 
I and Leo II) reveal that they too are closely aligned with 
J mean , then our results suggest that such a setup would be 
extremely unlikely in our ACDM galaxy formation model. 

Fig. [9] shows the probability that the orbital angular 
momentum of a random N out of 11 satellites falls close to 
the short axis of the disc-of-satellites or J mcan . The mean 
orbital angular momentum of N satellites could still lie close 
to the axis (as in the MW) , even if the angular momentum of 
the individual satellites are widely distributed. In Fig. HOh ,. 
we show the distribution of the angle between the mean 
orbital angular momentum of the 6 best aligned satellites, 
J mcan , and the short axis of the disc-of-satellites, c sat . 

Fig. llOfa . shows that there is an overabundance of 
systems that fall in the highest bin, corresponding to 
an angle of 9 < 25 degrees. The fraction of systems 
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Figure 10. The distribution of the angle between the mean or- 
bital angular momentum, J^ ean , of the 6 best aligned satellites 
and the short axis of the disc- of- satellites, c sat . We show the dis- 
tribution for both the model satellites (a) in the upper panel, and 
for the associated randomized dark matter particles (b) in the 
lower panel. The dashed red line indicates a uniform distribution 
in cos(6>). 



whose Jmcan li es within 25 degrees of 6 sat is around 17%, 
well above the 5% expected for a uniform distribution. 
This is not an exceedingly l arge p roportion, and although 
iMetz. Kroupa. fc Libeskindl (|2008h failed to find a single ex- 
ample in 6 simulations, our result - that such an alignment 
occurs roughly 17% of the time - is consistent with their 
findings. There is roughly a 50% chance of seeing such an 
alignment in just 1 out of 6 systems drawn at random from 
our sample. Yet, the non-uniformity of the distribution of 
the angle between J^ em and c sa t is statistically very signif- 
icant: a Kolmogorov- Smirnoff test rules out the hypothesis 
that this sample is drawn from a random distribution at 
very high confidence. 

We contrast these results with the histogram in the 
lower panel of Fig. 110b which shows the distribution of the 
same angle but now with Jm Can defined by the matched 
randomly selected control particles. In this case, the dis- 
tribution is consistent with uniformity: the KS probability 
of drawing such a sample from a uniform distribution is 
roughly 20%. In other words, the orbital angular momen- 
tum of the control particles (which were chosen to be at 
the same radial distances as the model galaxies), exhibits 



no preference for pointing near the normal of the plane they 
define. 



4 CONCLUSION AND DISCUSSION 

In this paper we have used a large cosmological TV-body 
simulation of the ACDM cosmology, populated with galax- 
ies using the Galform semi-analytic model of galaxy for- 
mation, to investigate the expected frequency of the pe- 
culiar alignments exhibited by the satelli te galaxies of the 
Milky Way. We c onfirm previous work fe.g | Kang et alj[2005l : 
iLibeskind et all 120051 : [ Zentner et al. 2005) which found 
highly flattened spatial distributions of satellites within the 
much more spherical dark matter distributions of simulated 
galactic halos. The satellite ellipsoid tends to be aligned with 
the shape of the dark matter halo, although in a small frac- 
tion of cases it is anti-aligned. 

We have also used our simulation in an attempt to un- 
derstand the dynamics of satellites within Milky Way type 
halos. We have found that systems that resemble the Milky 
Way - with at least 3 satellites orbiting within 30 degrees of 
the normal to the plane defined by the 11 brightest of them 
- occur roughly 35% of the time. Th i s is co nsistent with the 
work of lMetz. Kroupa. fc Libe skind ( 20081) who found 2 out 
of 6 galaxies showing this behavior. However, the fraction 
of systems which have 8 satellites all orbiting with angular 
momentum aligned to within 30 degrees of the normal to 
their plane is less then 1.5% at the 95% confidence level. 
Thus, if reliable proper motion estimates for the remaining 

5 "classical" Milky Way satellites reveal that they too orbit 
near this plane, this coincidence would be difficult to explain 
within our ACDM galaxy formation model. In roughly 17% 
of our simulated galactic halos we find systems in which the 
mean orbital angular momentum of the 6 best aligned satel- 
lites lies within 25 degrees of the short axis of the satellites' 
plane (and, in 30% of cases within 35 degrees), as seen in 
the Milky Way. 

The aligned satellites in the simulations do not make 
up a rotationally supported dis c-of-satellites. Similarly, and 
contra ry to the suggestion by IMetz. Kroupa. fc Libeskindl 
(2008), the rough alignment of some of the MW's satellites 
with the pole of their disc-of-satellites does not imply that 
the system is rotationally supported. This can only be ascer- 
tained by comparing the angular momenta of the satellites 
with the values for circular orbits of the same energy. Only 
this comparison (which requires various model assumptions) 
can expose how circular the orbits really are. If indeed the 
orbits of the Milky Way satellites turn out to be rotationally 
supported, this would be difficult to explain in the ACDM 
model. 

The main result of this paper is the robust quantifica- 
tion of the fraction of Milky Way type galaxies in the ACDM 
cosmology that resemble our Galaxy as far as its anisotropic 
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satellite distribution and dyna mics are concerned. While 



previous analysis of simulations ([Libeskind et al| 2005 2007; 
IZentner et ai1l2005l : iMetz. Kroupa. fc LibeskindlboOSI ) have 



come to some of the same conclusions as we do here, their 
statistical significance was compromised by small sample 
sizes (6, 3, 3 and 6 haloes respectively). Our large sample 
of over 400 Milky Way sized haloes allows us robustly to 
quantify the expected frequency of both satellite flattening 
and coherent satellite motion. 

Future surveys such as Pan-STA RRS or SkyMapper 
|Kaiser et all |2002| ; iKeller et ail I2007T I will obtain a more 
complete census of the Milky Way's satellite population, 
including measurements of the galaxy luminosity function, 
and better determinations of their kinematic properties and 
their relationship with the "great pancake" of bright s atel- 
lites . It was originally suggested bvlBenson et alj (|2002| ) (see 
also lKoposov et al.ll2007l ; lTollerud et al.ll2008l) . that new sur- 
veys may reveal hundreds, if not thousands, of faint small 
satellites. However, it seems highly unlikely that the known 
population of bright Milky Way satellites will change signif- 
icantly and thus that the conclusions of this paper - which 
pertain to the bright satellites - will be significantly altered. 
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